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Abstract 

The two particle interferometry method to determine the size of the emitting source after a heavy ion collision is extended. 
Following the extension of the method to spherical expansion dynamics, here we extend the method to rotating systems. 
It is shown that rotation of a cylindrically symmetric system leads to modifications, which can be perceived as spatial 
asymmetry by the ’’azimuthal HBT” method. 

We study an exact rotating and expanding solution of the fluid dynamical model of heavy ion reactions. We consider 
a source that is azimuthally symmetric in space around the axis of rotation, and discuss the features of the resulting 
two particle correlation function. This shows the azimuthal asymmetry arising from the rotation. We show that this 
asymmetry leads to results similar to those given by spatially asymmetric sources. 
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1. Introduction 

Two particle interferometry was adapted to heavy ion 
physics around 1980, and turned out to be a sensitive tool 
of determining the space-time extent of the source of the 
emission in heavy ion collisions. It was a significant early 
discovery that pion interferometry of exploding sources 
modifies the determination of the radius of the source [1]. 
This made the analysis a valuable tool for studying the 
dynamic expansion of nuclear sources. 

Heavy ion collisions with finite impact parameters create 
systems where we have a large net angular momentum in 
the initial state, which leads to a rotating [2 and expand¬ 
ing fireball. If the formed Quark-Gluon Plasma (QGP) 
has low viscosity [3], one can expect new phenomena like 
rotation or turbulence, which shows up in form of a start¬ 
ing Kelvin-Helmholtz instability (KHI) 0. Rotation in 
heavy ion collision recently has been considered, and here 
we study a new class of exact hydrodynamic solutions for 
three dimensional, rotating and expanding cylindrically 
symmetric fireball m- 

The created system in relativistic heavy-ion collisions is 
microscopic and short-lived so only the momentum spec¬ 
trum of the emitted particles can be measured directly. 
However, the space-time structure of the collision region 
can be studied using Hanbury-Brown-Twiss interferome¬ 
try 0. This technique uses two particle correlations p] , to 
probe the space-time shape of the particle emission zone. 
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The size and shape of reaction zone become thus accessi- 
ble, with the ’’azimuthal HBT” method |91fl5). 

The Differential Hanbury Brown and Twiss (DHBT) 
method has been introduced earlier in [8]. Previously the 
method has been applied to high resolution Particle in Cell 
Relativistic (PICR) fluid dynamical model [TO] results, and 
it was shown that rotation can be detected by this mod¬ 
ified method. The same PICR model was also used to 
calculate the vorticity of the flow, and due to the equipar- 
tition between the spin and orbit rotation the polarization 
of emitted particles was evaluated and turned out to 
be significant. This early prediction was verified recently 
|18j , by the experimental study of A and A polarization in 
Au+Au reactions in the energy range of ^/snn = 7.7—39 
GeV/nucl. Furthermore the A and A polarizations pointed 
in the same direction that verified the mechanical, equipar- 
tition origin of the polarization in contrast to electromag¬ 
netic origin, which would have led to opposite polarizations 
for A and A. 

In this work we calculate two pion correlation function 
for a rotating and expanding QGP, formed in Pb+Pb col¬ 
lisions at y/sNN = 2.76 TeV/nucl, and impact parameter 
b = 0.7 b max , by using the exact hydro model m we de¬ 
termine the effect of rotation on the correlation function 
(CF) for detectors at different positions. Finally we fit re¬ 
sults by ’’azimuthal HBT” to extract the apparent size of 
the rotating system in different directions. 

2. Correlation Function 

We use a simple Exact Model 7>] for expanding and ro¬ 
tating systems to demonstrate the sensitivity of the two 
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particle correlation method to diagnose rotation. Both 
polarization [T5] and the two particle correlation [50] were 
already evaluated for this model, with parametrizations 
adapted for peripheral heavy ion collisions [J>|. 

We consider an azimuthally symmetric system around 
the rotation axis, y, with Gaussian density profiles with 
characteristic radii, R and Y and constant temperature T. 
The initial parameters are given in Table [l] and the initial 
temperature is taken to be T = 200 MeV. 

The source function, S(x, k), giving the emission rate in 
the phase space, x, k , should be integrated over all points, 
x , of the emitting source to obtain the correlation function: 


t 

(fm/c) 

Y 

(fm) 

Y 

(c) 

UJ 

(c/fm) 

R 

(fm) 

R 

(c) 

(Rad) 

0 . 

4.000 

0.300 

0.150 

2.500 

0.250 

0.000 

3. 

5.258 

0.503 

0.059 

3.970 

0.646 

0.307 

8. 

8.049 

0.591 

0.016 

7.629 

0.779 

0.467 


Table 1: Time dependence of characteristic parameters of the fluid 
dynamical calculation presented in Ref. [6]. R is the average trans¬ 
verse radius, Y is the longitudinal length of the participant system, 
tp is the angle of the rotation of the interior region of the system, 
around the y-axis, measured from the horizontal, beam, z-direction 
in the reaction, [x, z], plane, R and Y are the speeds of expansion in 
transverse and longitudinal directions, and co is the angular velocity 
of the internal region of the matter. 
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Here the spatial integral is performed in cylindrical co¬ 
ordinates, s y , s p , ip, where s y , and s p are scaling variables, 
s v = V 2 1Y 2 and s p = (x 2 + z 2 )/R?. 

The correlation function was evaluated the same way 
as in Ref. [20] , We should see that the source function 
explicitly depends on the velocity field, v B , which includes 
both the expansion and the rotation of the system. 

As the rotation axis is the y-axis, the [x, z] plane is the 
reaction plane. Due to the azimuthal symmetry the radius 
of the system in the [x, z] plane is R. 

According to the conventions of two particle correlation 
functions in heavy ion physics, the z-axis is the beam axis 
and determines the LONG direction. In the transverse 
plane, the axaxis (the direction of impact parameter) is 
transverse to the beam direction. In this way the OUT 
direction is the ^-direction. The remaining y-axis deter¬ 
mines the SIDE direction. Due to the azimuthal symmetry 
of our specific model the results for the LONG and OUT 
directions should be identical. 

The velocity field in x, y, z coordinates is given by 

v s = [Ry/s^ sin(</>) + Ruy/s^ cos (tp), 

• _ • \ ( 2 ) 
Yy/s y ,Ry/s~ p cos {if) - RiOyfsp sin(<^)J , 

where w is the angular velocity, and ip is the angle of 
rotation. 

The mean transverse radius is R = \JXZ, and we use 
this value when the exact model is studied. 

In the practical calculations we use a detectors placed 
at k + = {k x ,k y ,k z )k = (0.924,0,0.383)fc, k = (1,0,0)fc 
and k~ = (0.924,0, —0.383)fc, which are orthogonal to the 
rotation axis, y. 


3. Results 

We calculate the CF for different values of the angular 
velocity, uj, to see how it is affected. The CFs are shown 
in Figs. [T] and [2| 
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Figure 1: (color online) Correlation Function, C(k,q), for the exact 
hydro model for the q = qouT direction. R = 2.50 fm, R = 0.25 c, Y 
= 4.00 fm, Y = 0.30 fm, uj = 0.30 c/fm, at t = 0.0 fm/c with k = 5 
fm -1 . The solid black line is for measuring the correlation function 
at k~ = (0.924, 0, —0.383)/c, the dashed red line is for k = (1,0, 0)fc 
and the dotted blue line is for k + = (0.924, 0, 0.383)Au 


Subsequently the correlation functions is fitted by the 
azimuthal HBT method and parametrization: 


C(q,k) = 1 +expi- ^ (MjR^ik) , (3) 

V i,j=L,o,s ) 

The obtained radius parameters in the directions L and 
O are identical due to the symmetry of the model, thus 
Roo = Rll- We took different fc-values. 

We can compare the different radius parameters, R 2 j, 
obtained by fitting the results of the CF obtained from 
the rotating and azimuthally symmetric system, to the 
’’azimuthal HBT” parametrization of eq. ([3|. The ob¬ 
tained values, R = Roo = Rll , are shown in Fig. [3j 

The CF increases for larger values of u> which corre¬ 
sponds to a decrease in the measured radius. The approx¬ 
imate radius decrease for to = 0 to 0.15 c/fm and oj = 0 to 
0.30 c/fm is 3-4% and 15% respectively. 

In Fig. [3] we can see how the measured radius depends 
on the rotation velocity u> at time t = 0 fm/c.. For a 
later times we would see a similar effect. Even though 
u> becomes smaller for later times, the expansion velocity 
and size of the system are both influencing the CF together 
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Figure 2: (color online) Correlation Function, C(k,q), for the exact 
hydro model, for the q = qouT direction. R = 2.50 fm, R = 0.25 
c, Y = 4.00 fm, Y = 0.30 fm at t = 0.0 fm/c with k = 5 fm -1 . 
The solid black line is for uj = 0.30 c/fm, the dashed red line is for 
L) = 0.15 c/fm and the dotted blue line is for co = 0.00 c/fm. 
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Figure 3: (color online) Ratio of radius from the fit for the correla¬ 
tion function in Fig. ^and[2] for different directions as a function of 
a;, the black line is for k~ = (0.924,0, —0.383)fc, the red line is for 
k = (1,0, 0)fc and the blue line is for k+ = (0.924, 0, 0.383)fc. Rq is 
the observed radius of the system without rotation. 

with the rotation. On the other hand there is no effect on 
the radius parameters if either the rotation or expansion 
is zero [5U]j . 

A higher rotation velocity will decrease the measured 
size of the system, it also decreases more rapidly for larger 
values of ui as can be seen from the slope going from u> = 
0 to 0.15 c/fm and oj = 0.15 to 0.30 c/fm in Fig. [3] 
Asymmetry in the size is present if measured at different 
directions if the system is rotating. If the rotation were 
reversed the correlation function will also change, where 
the black and blue lines in Figs. [T] and [3] are exchanged. 

The detector at k + = (0.924, 0, 0.383)A: shows a smaller 
measured radius for the exact hydro model while the ra¬ 
dius showed at k~ = (0.924,0, —0.383)fc is larger. This 
is also dependent on expansion velocity, temperature and 
size of the system. The axial size, Y , is not affected by the 
rotation. 

Thus the model results show that rotation significantly 


influences the HBT evaluation similarly like the expansion, 
which influences data significantly, e.g. in publications 

mm- 


4. Conclusion 

Different values of the angular velocity will change the 
measured ” azimuthal HBT” size parameters of the system. 
It will also create smaller and larger values for the cor¬ 
relation function when measuring at different directions. 
That the cylindrically symmetric system is rotating will 
be observed as an asymmetric object in ’’azimuthal HBT” 
analysis. 


Acknowledgements 

This work is supported by the Research Council of Nor¬ 
way, Grant no. 231469. 

References 

[1] Scott Pratt, Phys. Rev. Lett. 53, 1219 (1984). 

[2] L.P. Csernai, V.K. Magas, H. Stocker, and D.D. Strottman, 
Phys. Rev. C 8, 024914 (2011). 

[3] L.P. Csernai, J. I. Kapusta, L. D. McLerran, Phys. Rev. Lett. 
97, 152303 (2006). 

[4] L.P. Csernai, D.D. Strottman, and Cs. Anderlik, Phys. Rev. C 
85, 054901 (2012). 

[5] T. Csorgo and M.I. Nagy, Phys. Rev. C 89, 044901 (2014). 

[6] L. P. Csernai, D. J. Wang, and T. Csorgo, Phys. Rev. C 90, 
024901 (2014). 

[7] M. Nagy, V. Konye, and T. Csorgo, Talk at the Balaton Work¬ 
shop, July 11-17, 2015, Tihany, Hungary. 

[8] L.P. Csernai and S. Velle, Int. J. Mod. Phys. E 23, 1450043 
(2014). 

[9] S. A. Voloshin and W. E. Cleland, Phys. Rev. C 53, 896 (1996); 
Phys. Rev. C 54, 3212 (1996). 

[10] U. A. Wiedemann, Phys. Rev. C 57, 266 (1998). 

[11] M. A. Lisa, U. Heinz and U. A. Wiedemann, Phys. Lett. B 489, 
287 (2000) 

[12] U. Heinz and P. F. Kolb, Phys. Lett. B 542, 216 (2002). 

[13] U. Heinz, A. Hummel, M.A. Lisa and U.A. Wiedemann Phys. 
Rev. C 66, 044903 (2002). 

[14] F. Retiere and M. A. Lisa, Phys. Rev. C 70, 044907 (2004). 

[15] G. Graef, M. Bleicher and M. Lisa Phys. Rev. C 89, 014903 
(2014). 

[16] L.P. Csernai, S. Velle, and D. J. Wang, Phys. Rev. C 89, 034916 
(2014). 

[17] F. Becattini, L. P. Csernai, and D. J. Wang, Phys. Rev. C 88, 
034905 (2013). 

[18] M. A. Lisa, Invited talk at the XI Workshop on Particle Cor¬ 
relations and Femtoscopy (WPCF2015), November 3-7, 2015, 
Warsaw, Poland. 

[19] Yilong Xie, Robert C. Glastad, and Laszlo P. Csernai, Phys. 
Rev. C 92, 064901 (2015) 

[20] S. Velle and L.P. Csernai, Phys. Rev. C 92, 024905 (2015). 

[21] L.Q. Shan, F.J. Wu, J.L. Liu, et al., J. Phys. G 36, 055003 
(2009). 

[22] Q.F. Li, M. Bleicher, and H. Stocker, J. Phys. G 34, 2037 (2007). 

[23] S. Pratt, and J. Vredevoogd, Phys. Rev. C 78, 054906 (2008). 

[24] W.N. Zhang, Y.Y. Ren, C.Y. Wong, Phys. Rev. C 74, 024908 
(2006). 

[25] D.K. Srivastava, Pliys. Rev. C 71, 034905 (2005). 


3 










